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A long-standing aim of quantum information research is to understand what gives quantum com¬ 
puters their advantage. This requires separating problems that need genuinely quantum resources 
from those for which classical resources are enough. Two examples of quantum speed-up are the 
Deutsch-Jozsa and Simon’s problem, both efficiently solvable on a quantum Turing machine, and 
both believed to lack efficient classical solutions. Here we present a framework that can simulate 
both quantum algorithms efficiently, solving the Deutsch-Jozsa problem with probability 1 using 
only one oracle query, and Simon’s problem using linearly many oracle queries, just as expected of 
an ideal quantum computer. The presented simulation framework is in turn efficiently simulatable 
in a classical probabilistic Turing machine. This shows that the Deutsch-Jozsa and Simon’s problem 
do not require any genuinely quantum resources, and that the quantum algorithms show no speed-up 
when compared with their corresponding classical simulation. Finally, this gives insight into what 
properties are needed in the two algorithms, and calls for further study of oracle separation between 
quantum and classical computation. 


Quantum computational speed-up has motivated much 
research to build quantum computers, to find new algo¬ 
rithms, to quantify the speed-up, and to separate clas¬ 
sical from quantum computation. One important goal 
is to understand the reason for quantum computational 
speed-up; to understand what resources are needed to do 
quantum computation. Some candidates for such neces¬ 
sary resources include superposition and interference [1], 
entanglement [2], nonlocality [3], contextuality [4-6], and 
the continuity of state-space [7]. In this paper we look at 
two so-called oracle problems: the Deutsch-Jozsa prob¬ 
lem [8, 9] and Simon’s problem [10, 11], thought to show 
that quantum computation is more powerful than clas¬ 
sical computation. There are quantum algorithms that 
solve these problems efficiently, and here we present a 
framework that can simulate these quantum algorithms. 
Quantum Simulation Logic (QSL), that in itself can be ef¬ 
ficiently simulated on a classical probabilistic Turing ma¬ 
chine. To underline this, we provide Python implemen¬ 
tations of these simulations as supplementary material, 
that give the correct output in seconds for 100,000-qubit 
inputs to the Deutsch-Jozsa algorithm and 512-qubit in¬ 
puts to Simon’s algorithm, for random choice of the func¬ 
tion under consideration. This shows that no genuinely 
quantum resources are needed to solve these problems 
efficiently, but also tells us which properties are actually 
needed: the possibility to choose between two aspects 
of the same system within which to store, process, and 
retrieve information. 

The quantum algorithm for the Deutsch-Josza problem 
has been used extensively for illustrating experimental 
realizations of a quantum computer, while Simon’s algo¬ 
rithm served as inspiration for the later Shor’s algorithm 
[13]. Both problems involve finding certain properties of 
a function /, and this can be done efficiently in a quan¬ 
tum computer given a quantum gate that implements 
the function, known as an oracle. When using a classi¬ 
cal function as an oracle (an oracle that is only allowed 
to act on classical bits as input and giving classical bits 



FIG. 1: Circuits for the Deutsch-Jozsa algorithm. A. The 
algorithm itself. This circuit uses an n-qubit query-register 
prepared in the state jO), and a target prepared in jl). It pro¬ 
ceeds to apply Hadamard transformations to each qubit. The 
function / is embedded in an oracle Uf, and this is followed 
by another Hadamard transformation on each query-register 
qubit. The measurement at the end will test positive for jO) 
if / was constant, and negative if / was balanced. B. Oracle 
construction for balanced functions. The arbitrary permuta¬ 
tion TT = Itt)®)) { x \ of the computational basis states is 

constructed from CNOT and Toffoli gates [12] (rr”^ = ttQ. 
At the center is a CNOT gate from the most significant qubit 
\xn) to the target \y). With rr as the identity permutation, 
the oracle performs the balanced function f' that is 1 for all 
inputs with the most signihcant bit set. Any other balanced 
function f{x) = f'{TT{x)) can now be generated by choosing 
a different tt. For a function that is constant zero, the CNOT 
gate is omitted, while the constant one function should re¬ 
place the CNOT with a Pauli-X gate acting on the target. 


as output), the Deutsch-Jozsa problem has an efficient 
solution on a classical probabilistic Turing machine [8, 
9], but for Simon’s problem it has been proven [10, 11] 
that no efficient solntion exists. In both the quantum 
and classical case the oracle is treated as a black box, 
i.e., its internal structure is not visible from the outside, 
and each oracle call counts as one operation; one unit 
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FIG. 2: QSL circuitry. A. QSL construction of a CNOT gate 
(computational bit in black, phase bit in grey). B. CNOT 
gate identity valid both for qubit and QSL gates. C. QSL 
construction of a Toffoli gate. D. QSL construction of a gen¬ 
eral permutation. 


of time. An algorithm that requires at most a polyno¬ 
mial number of oracle calls and at most a polynomial 
amount of additional resources is regarded as efficient, 
but the oracle itself does not need to be efficiently com¬ 
putable. In essence, the question asked is “If this oracle 
was efficiently computable, would the whole algorithm be 
efficiently computable?” For both the mentioned prob¬ 
lems, the answer is “yes” in quantum computation. For 
Simon’s problem, the answer is “no” for the classical- 
function oracle, giving evidence for a separation relative 
to the oracle between quantum computation and classi¬ 
cal computation. For the Deutsch-Jozsa problem such a 
separation is only obtained by requiring a deterministic 
solution. 


I. THE DEUTSCH-JOZSA ALGORITHM 
AND QUANTUM SIMULATION LOGIC 

The Deutsch-Jozsa problem is the following: Suppose 
that you are given a Boolean function f{x) : {0,1}” >-)■ 
{0,1} with the promise that it is either constant or bal¬ 
anced. The function is constant if it gives the same out¬ 
put (1 or 0) for all possible inputs, and it is balanced 
if it gives the output 0 for half of the possible inputs, 
and 1 for the other half. Your task is now to distin¬ 
guish between these two cases [9, 14]. Given such a func¬ 


tion, a classical Turing machine can solve this problem 
by checking the output for 2"“^ + 1 values of the in¬ 
put; if all are the same, the function is constant, and 
otherwise balanced. A stochastic algorithm with k ran¬ 
domized function queries gives a bounded error proba¬ 
bility [9] less than 2^“^, showing that the problem is in 
the complexity class BPP (Bounded-error Probabilistic 
Polynomial-time solvable problems). 

A quantum computer obtains the function as a unitary 
that implements the function, a quantum oracle. For 
later convenience we use the gate implementation of the 
oracle Uf shown in Figure lb. This device will give the 
expected outcomes for regular function queries, i.e., in¬ 
vocations that reveal a single function value. The desired 
input number k should be inserted in the query-register 
as the quantum state \k) in the computational basis, and 
the answer-register should be in a computational basis 
state (say |0)). Applying the oracle and measuring the 
answer-register will reveal the output of /: if the answer- 
register was flipped, then the function value is 1 for that 
input k. Given this quantum oracle the Deutsch-Jozsa 
algorithm [9, 14] can solve the problem with a single 
query (see Figure la). In this algorithm Hadamards are 
applied to the query-register state, creating a quantum 
superposition of all possible inputs to the oracle. When 
the oracle is applied, this superposition is unchanged if 
the function is constant, and changed to a different su¬ 
perposition if the function is balanced. A change can be 
detected by again applying Hadamards and measuring 
in the computational basis. In the ideal case the error 
probability is zero, showing that the problem is in EQP 
(Exact or Error-free Quantum Polynomial-time solvable 
problems [15]). 

The QSL versions of the quantum algorithm will use 
the same set-up as in Figure 1, but instead of qubits 
it will use pairs {bz,bx) containing two classical bits, a 
“computational” bit b^ and a “phase” bit bx- These con¬ 
stitute the elementary systems of our QSL framework, 
the “QSL bits”. State preparation of a computational 
single qubit state \k) is associated with preparation of 
{k, X) where A is a random evenly distributed bit. Mea¬ 
surement in the computational basis is associated with 
readout of the computational bit followed by randomiza¬ 
tion of the phase bit, somewhat like the uncertainty rela¬ 
tion or really measurement disturbance as seen in quan¬ 
tum mechanics. Accordingly, measurement of the phase 
bit will be followed by a randomization of the computa¬ 
tional bit. These constructions of state preparation and 
measurement prohibits exact preparation and readout of 
the system; the upper limit is one bit of information per 
QSL bit {bz, bx). 

It is relatively simple to simulate the single qubit gates 
used in the quantum algorithms: an X gate inverts the 
computational bit (preserving the value of the phase bit 
bx), a, Z gate inverts the phase bit (preserving bz), and a 
H gate switches the computational and phase bits. These 
are constructed so that the resulting gates obey the quan¬ 
tum identities XX = ZZ = HH = / and HZH = X. It 
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is also possible to define QSL V and “phase” gates, but 
we will not discuss these and their associated identities 
here as they are not needed for the present task. A sim¬ 
ple circuit is to prepare the state |0), apply an ff gate, 
and measure in the computational basis. This will give 
a random evenly distributed bit as output. 

A system of several qubits can now be associated with 
a system of several QSL bits containing two internal bits 
each. A QSL simulation of the quantum Controlled- 
NOT (CNOT) gate can be constructed from two clas¬ 
sical reversible-logic CNOTs. One classical CNOT con¬ 
nects the computational bits in the “forward” direction, 
while the other connects the phase bits in the “reverse” 
direction, see Figure 2a. This again is constructed to 
enable use of the same identities as apply for the quan¬ 
tum CNOT, most importantly a phenomenon sometimes 
called “phase kick-back” in quantum computation [14], 
see Figure 2b. 

Our QSL framework is inspired by a construction of 
Spekkens [16] that captures many, but not all, properties 
of quantum mechanics. So far, the presented framework 
is completely equivalent to Spekkens’ model (the H gate 
appears in a paper by Pusey [17]). Both frameworks are 
efficiently simulatable on a classical Turing machine (see 
above and Ref. [16]). Furthermore, both have a close re¬ 
lation to stabilizer states and Clifford group operations 
[17], and perhaps more interestingly, both frameworks 
capture some properties of entanglement, enabling pro¬ 
tocols like super-dense coding and quantum-like telepor¬ 
tation, but cannot give all consequences of entanglement, 
most importantly, cannot give a Bell inequality violation. 

Returning to the Deutsch-Josza problem, for n = 1 the 
oracle for the balanced function f{x) = x is a CNOT, 
which means that Figure 2b contains the whole Deutsch- 
Josza algorithm when using the prescribed initial state 
|0)|1). With this initial state, the query-register mea¬ 
surement result would be 1, so that this oracle gives the 
same output from the QSL algorithm as from the quan¬ 
tum algorithm. Indeed, all oracles for n = 1 and n = 2 
only use CNOT and X gates, so their simulation will 
give the same behaviour as the quantum oracles. For 
n = 1 and n = 2, the Deutsch-Jozsa algorithm is already 
known to have an implementation that does not rely on 
quantum resources [18], and also, the gates used so far 
can be found within the stabilizer formalism, and here 
the Gottesmann-Knill theorem tells us that they can be 
simulated efficiently [19]. 

However, for n > 3, the Deutsch-Jozsa oracle needs 
the Toffoli gate, and since the Toffoli gate cannot be 
efficiently simulated using the stabilizer formalism [19], 
nor is present in Spekkens’ framework [16, 20], it has so 
far been believed that the Deutsch-Jozsa algorithm does 
not have an efficient classical simulation. Here, we need 
to point out that our task is not to create exact Tof¬ 
foli gate equivalents, or even simulate the full quantum- 
mechanical system as such. It suffices to give a working 
efficient QSL version of the Deutsch-Jozsa algorithm. We 
therefore choose not to represent Toffoli gates exactly. 


but design the gate so that it implements a classical Tof¬ 
foli in the computational bits, and some other gate in the 
phase bits. For simplicity we choose the identity map for 
the phase bits, see Figure 2c. 

By representing a general computational-state permu¬ 
tation TT in the same manner as the Toffoli (see Fig¬ 
ure 2d), the quantum oracle can immediately be trans¬ 
lated into an oracle within the QSL framework. The 
balanced-function oracle in Figure lb leaves the phase 
bits unchanged except for the one belonging to the most 
significant QSL bit in the query-register. A constant- 
function oracle leaves all phase bits unchanged. Since 
Hadamard gates are included before and after the ora¬ 
cle, measurement of the computational bits will reveal 
the oracle’s effect on the phase bits of the query-register, 
solving the problem by only one oracle query. 

This QSL algorithm uses equally many QSL bits and 
QSL gates as the quantum algorithm uses qubits and 
quantum gates. The time and space complexity are 
therefore identical to the complexity of the quantum algo¬ 
rithm. The QSL algorithm can be efficiently simulated on 
a classical probabilistic Turing machine, using two clas¬ 
sical bits for each simulated qubit and at most twice as 
many classical reversible gates as quantum gates. The er¬ 
ror probability is 0, and in fact, the same correct output 
is generated for all possible random values in the proba¬ 
bilistic machine, the machine’s random value could there¬ 
fore be replaced with a constant value without destroying 
the simulation. Therefore, this QSL algorithm can be 
simulated on a classical deterministic Turing machine, 
showing that, relative to the oracle, the Deutsch-Jozsa 
problem is in P (Polynomial time solvable problems). 


II. SIMON’S PROBLEM 

Simon’s problem is to decide whether a function /(x) : 
{0,1}" I—)■ {0,1}" is one-to-one or two-to-one invariant 
under the bitwise exclusive-OR operation with a secret 
string s (or equivalent, invariant under a nontrivial XOR- 
mask s [10, 11]). In both cases, /(x) = /(x') if and only 
if x' = X0S (bit-wise addition modulo 2); if s = 0 then / 
is one-to-one and otherwise / is two-to-one. When using 
a classical function as an oracle, a classical probabilistic 
Turing machine needs an exponential number of oracle 
evaluations for this task [10, 11], showing that relative to 
this oracle, Simon’s problem is not in BPP. 

Simon’s quantum algorithm can distinguish these two 
cases in an expected linear number of oracle queries [10, 
11], showing that the problem is in BQP (Bounded-error 
Quantum Polynomial time [15]). A circuit diagram of 
the quantum subroutine and one realization of the oracle 
used in Simon’s algorithm is shown in Figure 3a-b. An 
explicit construction of the gate Ug can be obtained by 
choosing a basis {t'*} for the part of the binary vector 
space orthogonal to s. If s = 0, the basis consists of n 
vectors, otherwise n — 1 vectors. For every entry vj’ = 1, 
connect a CNOT from query-register bit j to ancilla- 
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FIG. 3: Circuits for the quantum subroutine of Simon’s al¬ 
gorithm. A. The subroutine itself. The n-qubit input and 
output registers are prepared in the state |0)|0). Apply 
Hadamard transformations to the query-register, the function 
/ embedded in an oracle Uf, and finally another Hadamard. 
The measurement at the end will give a random bit sequence y 
such that y • s = 0 (mod 2). B. Oracle construction. The ora¬ 
cle Uf uses an internal n-qubit ancilla, and the CNOT denotes 
a string of bit-wise CNOTs from each ancilla bit to the corre¬ 
sponding answer-register bit. The unitary Ua is constructed 
from CNOTs [21] (see the main text) so that it implements 
one representative function for each specific secret string s, 
and the permutation tt gives access to all functions /. The 
dotted gate combination Uf implements |a:)|0) i—>■ |a:)|/( 2 ;)) as 
required in Simon [10, 11], and the additional circuitry in Uf 
is there to map |2;)|2) to \x)\z (B f{x)) and reset the ancilla 
[22]. C. The modified subroutine for the deterministic algo¬ 
rithm. The n-qubit query- and answer-register are prepared 
in the state |0)]ui). The circuit applies Hadamard transfor¬ 
mations to the query- and answer-registers, the function / 
embedded in an oracle Uf, and finally another Hadamard be¬ 
fore measurement. 

register k (in realizations, it is beneficial to choose to 
minimize the number of CNOTs [21]). The permutation 
after Us enables all possible functions. 

Let us illustrate the construction of Ug with the fol¬ 
lowing example. Suppose s = (1,0,1), then we choose 
the basis = {(1,0,1), (0,1,0)} which spans the bi¬ 
nary vector space orthogonal to s (note that the scalar 
product is defined mod 2). Now we see that t/(i,o,i) can 
be constructed from three CNOTs, one from the first 
query-register qubit to the first ancilla qubit, one from 
the third query-register qubit to the first ancilla qubit, 
and one from the second query-register qubit to the sec¬ 
ond ancilla qubit. 

The quantum subroutine outpnt is a random bit vec¬ 
tor y such that y ■ s = 0 mod 2 (bit-wise dot product), 
uniformly distributed over all such bit vectors y. After 
an expected number of iterations that is linear in n, the 
subroutine will have produced n — 1 linearly independent 
values of y [10, 11], so that the resulting linear system of 
equations for s can be solved, giving one non-trivial solu¬ 
tion s*. If the function is one-to-one, this solution is just 
a random bit sequence, but if the function is two-to-one 


the solution is the actual secret string s. Querying the 
function for /(O) and /(s*) will give equal values if the 
function is two-to-one, and unequal values if the function 
is one-to-one, solving Simon’s problem. Note that as a 
consequence of solving Simon’s problem one also obtains 
the secret string s. 

The QSL versions of the subroutine and oracle is 
again obtained by substituting the quantum gates with 
their corresponding QSL gates. The QSL CNOT and 
Us constructions give the explicit map {x,p),(0,a) !—>■ 
{x,p® g'{a)),{f'{x),a) where 

fj{x) = X ■ (mod 2) = fj{x 0 s), 
s'M) = (mod 2). 

k 

The permutations tt and do not influence the phase, 
so the complete oracle Uf is the map (x,p),{z,w) !—>■ 
{x,p® g{w)),{^z® f{x),w) where 

f{x)=TT{f{x)) = f{x®s), 
g{w) = g'{a) © g'{w © a) = g'{w) 

By the use of Hadamard gates, the simulation of Simon’s 
subroutine sets the query-register phase entering the or¬ 
acle to p = 0”, while the answer-register phase w will 
be uniformly distributed. Measurement of the computa¬ 
tional bits after the final Hadamard will therefore give 
a random bit vector y = g{w) that is orthogonal to s, 
uniformly distributed over the possible values. This re¬ 
produces the quantum predictions exactly, showing that 
the simulation will work with this subroutine to solve Si¬ 
mon’s problem with the same expected linear number of 
iterations in n. Again, efficient simulation of the QSL 
framework on a classical probabilistic Turing machine 
shows that, relative to the oracle, Simon’s problem is 
in BPP. 

The derandomized quantum algorithm for Simon’s 
problem [23, 24] (with zero error probability), that shows 
that Simon’s problem relative to the quantum oracle is 
in EQP, is not immediately usable because it incorpo¬ 
rates a modified Grover search [25] which is not known to 
have an implementation within the QSL framework. The 
quantum algorithm is intended to remove already seen y 
and, crucially, prevent measuring ?/ = 0 which would be 
a failed iteration of the algorithm. The modified Grover 
search is presented [23] as applied to the output oi Uf, 
but since the final step of the Grover search is again Uj, 
the whole combination can be viewed as Uf preceded by 
another quantum algorithm. This other quantum algo¬ 
rithm, in a sense, finds an input to Uf such that the 
coefficient of the term |0) of the output is 0. 

Since we do not have a Grover equivalent in our frame¬ 
work, we would like a simpler solution. One way to 
achieve this is to prepare the initial answer-register state 
as the Hadamard transformation of a chosen state \w), 
see Figure 3c. Iterating \ w) over a basis for the bit-vector 
space will produce n values of y that span the bit vector 
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space orthogonal to s. Either this spans the whole space 
(s = 0), or gives a linear system of equations that has one 
non-trivial solution, equal to s. This QSL algorithm gives 
deterministic correct outputs just as the QSL Deutsch- 
Jozsa algorithm, and can therefore also be efficiently sim¬ 
ulated on a classical deterministic Turing machine, show¬ 
ing that relative to the oracle, Simon’s problem actually 
lies within P. 

This linear search procedure works in our setting be¬ 
cause the QSL permutation gate tt does not influence 
the phase bits. In the quantum-mechanical case, the 
quantum gate implementing tt does change what is some¬ 
times called phase information, e.g., in the simple case 
when the quantum tt is a CNOT. A more convoluted 
example that does not correspond to a simple modifi¬ 
cation of the phase information is when the quantum tt 
is a (non-stabilizer) Toffoli gate. Thus, the derandom- 
ized quantum algorithm [23] requires the more advanced 
modified Grover’s search algorithm. However, note that 
if the oracle is the smaller gate combination Uj in Fig¬ 
ure 3b that obeys only the original [10, 11] requirement 
ja;)j0) I—)■ ja;)j/(a;)) (for z = 0), then the quantum ver¬ 
sion of the simpler algorithm proposed above will work 
as well, avoiding the modified Grover search. 


III. CONCLUSIONS 

In conclusion, we have devised QSL equivalents of the 
Deutsch-Jozsa and Simon’s quantum algorithms. In the 
quantum algorithm, you are given a quantum oracle, a 
unitary gate that implements / by acting on qubits, your 
task is to distinguish two or more families of functions. In 
the presented QSL algorithms, you are given a gate that 
implements / by acting on QSL systems, and the same 
task. Using the QSL framework, we obtain the same suc¬ 
cess probability and the same time and space complexity 
as for the quantum algorithms, and these QSL algorithms 
are in turn efficiently simulatable in classical Turing ma¬ 
chines (implementations that can be run with very large 
inputs are provided as Supplementary material). There¬ 
fore, the Deutsch-Jozsa and Simon’s algorithms cannot 
anymore be used as evidence for a relativized oracle sep¬ 
aration between EQP and BPP, or even P. Both prob¬ 
lems are, relative to the oracles, in fact in P. 

Apparently, the resource that gives these quantum al¬ 
gorithms their power is also available in QSL, but not 
when using an oracle in the form of an ordinary classical 
function. Each qubit is simulated by two classical bits in 
a classical probabilistic Turing machine, and in the above 
oracles one of these bits is used for computing the func¬ 
tion, while the other is changed systematically when Uf 
is applied. This change (in phase information) can then 
be revealed by choosing to insert Hadamards before and 
after the oracle. This choice is present in both quantum 
systems and our simulation, and enables revealing either 
function output or function structure, as desired. Our 
conclusion is that the needed resource is the possibility 


to choose between two aspects of the same system within 
which to store, process, and retrieve information. 

While we still believe that quantum computers are 
more powerful than classical computers, the question 
arises whether oracle separation really can distinguish 
quantum-computational complexity classes from classi¬ 
cal complexity classes. There are many other examples 
of quantum-classical oracle separation, [26-28] and some 
simple cases can be conjectured to have efficient simu¬ 
lations in the QSL framework (e.g., 3-level systems [16] 
for the three-valued Deutsch-Jozsa problem [26]), but in 
general the question needs further study. It is possi¬ 
ble that the same technique can be used for direct com¬ 
putation using other quantum algorithms [13, 25], but 
this will take us out of the oracle paradigm. Also, gen¬ 
eral quantum computation has been conjectured to use 
genuinely quantum properties such as the continuum of 
quantum states, or contextuality [5, 6], which both are 
missing from the QSL framework [16]. In any case, nei¬ 
ther the Deutsch-Jozsa nor Simon’s algorithm needs gen¬ 
uinely quantum mechanical resources. 
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Appendix A: Python library for the Quantnm 
Simulation Logic 

#! /usrAiin/env p 3 rthon 

# Qiiantum_Simulation_Logic.py 

# Ve 2 rsion 1.0 

# Copyri^t Jan-Ake Larsson <jan-ake.larsson(91iu.s^ 2015 

inport random 
from math hiqxiirt log 

class QSL_syst€m(list): 
def_init_(self,value,n=None,phase=^one): 

# QSL states are either initialized fron bit strings or integers, and 

# preferrably the length n should be fixed if an integer is given 

list._init_(self, [0,0]) 

self .n^ 

if type (value) =^tr or type(value)=^inaryVector: 
self [0] ^nt (value, 2) 
if not self.n: 
self. n=len (value) 
else: 

self [0]=value 
if not self.n: 

self .n=int(log(valuet=1.0,2))+l # at least 1 bit 
self [l]^hase 
if not self[1]: 

self [Ugandan.randint(0, (l«self .n)-l) 

def_add__(self,other): 

new_canpF (self [0] «other. n) +other [0] 
new_phas€F (self [1] «other.n) bother [1] 
new_n=Belf .n-tother .n 

return QSL_system(new_conq5,new_n,new_phase) 

def_str__(self): 

coiTpf=4Din(self [0]) [2:] 
phase=^in(self [1] ) [2:] 

return "["+"0"*(self.n-len(canp))+conq>i-", "\ 

+"0" * (self. n-len(phase)) -tphaseh"]" 
def extend(self,other): 

self [0] = (self [0] «other. n) -hDther [0] 
self [1] = (self [1] «other. n) -hDther [ 1] 
self. n^elf. n+other. n 
def remove(self,n): 
self [0]/=l«n 
self [1] /=l«n 
self.n-=n 


def X(state,k): 

# k is the bit number to apply the gate on 
state[0] ^tate[0] ^ (l«k) 

def Z(state,k): 

# k is the bit number to apply the gate on 
state[l] ^tate[l] ^ (l«k) 

def Y(state,k): 

# k is the bit number to apply the gate on 
state[0] ^tate[0] ^ (l«k) 

state[l] ^tate[l] ^ (l«k) 
def H(state,k): 

# If equal do nothing, but if not equal, flip both 
if state[0]&(l«k) != state[l]&(l«k): 

Y(state,k) 
def CN(state,k): 

# k is the tuple (control bit number, target bit number) 

# Caiputational bit 

if state[0]&(l«k[0]): # If control ccarputational bit is 1 
state[0]=state[0]'‘(l«k[l]) # Flip target ccraputational bit 

# Phase bit 

if state[l]&(l«k[l]): # If target phase bit is 1 
state[l]=state[l]'‘(l«k[0]) # Flip control phase bit 

# Measurement output is rudimentary binary vectors 
class BinaryVector(str): 

# Lnmutable parent class: replace __new_, not_init_ 

def __new__(self ,value,n^one): 

if type (value) !^tr: 

retval=^ in (value) [2:] 
if nl^one: 

valu^"0"* (n-len (value) )+value 
return str.__new__(self,value) 
def_add__(self,other): 


retval^inaryVector(bin(int (self, 2) '‘int (other, 2)) [2: ], 
max(len(self) ,len(other))) 

return retval 

def_radd__(self,other): 

return self+other 
def dot(self,other): 

return bin(int(self,2)&Lnt(other,2)) .count("l")Z2 

def measureZ(state,l^one): 

# Measure in the ccmputational basis 
if l^^one: 

Ip^ange (St at e. n) 
retval="" 

for j in reversedOO: # Most significant bit first in string 
state[l]=state[l] ~ (random.randint(0, l)«j) 
if state[0]&(l«j): 

retval+="l" 

else: 

retval+="0" 

return BinaryVectorCretval) 

# Efficient representation of a random permutation _and_ its inverse, 
class RandPerm(dict): 

# This would be a "reversible random oracle" in conpiter science: if the 

# input has not been seen before, draw a randan output. If the input has 

# been seen, repeat the corresponding output. Make sure the random oracle 

# one-to-one, and provide the inverse, 

def_init_(self, n, parent?4Ione): 

diet._init_(self) 

self.N=2**n-1 
if parent^^one: 

self. inverse=flandPerm(n,self) 
else: 

self. inverse=parent 
def __getitem__(self, key): 

if (typeCkey) !=int and type(key) !=long) or key<0 or ke 3 S>self .N: 

raise ValueError 
try: 

retval = diet, g etitem (self, key) 
except: 

# Draw new items randomly, when needed, 
r et val^andon. randint (0, self. N) 
while retval in self .valuesO: 

ret val=random. randint (0, self. N) 
self [key] retval 
self. inverse [retval] %ey 
return retval 

def Peniiute(state,perm,k): 

# perm is the permutation of the number the computational bits form 

# k is a list or tuple of bit indices to apply the gate on 
x=0 

for i,j in enumerate(k): 
if state[0]&(l«j): 
x+=(l«i) 

# Bitwise xor to conpare input and output bits of the permutation 
y^sx^permW 

# Change the bits that should be changed 
for i,j in enumerate(k): 

if y & (l«i): 

X(state,j) 

def InversePermute(state,perm,k): 
try: 

inv=perm. inverse 
except: 

inv=[perm.index(i) for i in xrange(len(pernO)] 

Permute(state, inv,k) 



Appendix B: Python simulation of the QSL 
Deutsch-Jozsa algorithm 


Appendix C: Python simulation of the QSL Simon’s 
algorithm 


#! /usrAiin/env p 3 rthon 

# QSL_Deutsch_Jozsa.py 

# Version 1.0 

# Copyri^t Jan-Ake Larsson <jan-ake.larsson(91iu.s^ 2015 

from Quantum_Siiiiulation_Logic in^xirt ♦ 

def Ufbalanced(state,pernO: 

Peimute(state,perm,xrangeCl,state.n)) 

CN(state, (state.n-1,0)) 

InversePeimuteCstate,perm,xrangeCl,state, n)) 

def UfconstantO(state,perm): 

Peimute(state,perm,xrange(l,state.n)) 
pass 

InversePeimuteCstate,perm,xrange(l,state, n)) 

def Ufconstantl(state,perm): 

Peimute(state,perm,xrange(l,state.n)) 

X (state, 0) 

InversePeimuteCstate,perm,xrange(l,state, n)) 

### With n=20, printing all function values takes considerable time. 

### With n=100, the number of states cannot even be input in xrangeO . 

### Even with n=100000, the below Deutsch-Jozsa algorithm answers in seconds. 
n=5 

### Choose an oracle: one of UFbalanced, UfconstantO or Ufconstantl, and a 

### random permutation (can be explicit and non-randon if desired, see below) 

Uf=Ufbalanced 

^ftJf=UfconstantO 

^ftjf <Jf constant 1 

pen[i=RandPem(n) 

### An ejqjlicit peimutation can also be used, for example mimi cking a 
### 4-Toffoli as follows here. 

# pemF=[0,1,2,3,4,5,6,15,8,9,10,11,12,13,14,7] 

### This representation of the peimutation is internal to the oracle. Note 
### that the efficiency of the Deutsch-Jozsa algorithm refers to uses of the 
### oracle, not the oracle storage or internals. The imjjortant feature is that 
### the oracle is used only once. 

print n, "bit input" 

# For small n, the ordinary classical algorithm can be used 
if n<16: 

print "All function values" 
for i in xrange(2*’tn): 
state=QSL_syst€sn(i ,n)-H3SL_system(0,1) 

Uf (state,perm) 

outcc*ne=measureZ(state) 

print "f(Zs)=y,s"“/o(outcaiie[:-l] ,outcome[-l]) 

print "Deutsch-Jozsa:", 
state=QSL_syst€m(0 ,n)-H5SL_system( 1,1) 
for i in xrange(state.n): 

H(state,i) 

print "(one oracle call only)", 

Uf (state,pemO 

for i in xranged,state.n): 

H(state,i) 

outcc*ne^easureZ(state,xrange(l,state .n)) 
if outcom^="0"’tn: 

print "Constant" 
else: 

print "Balanced" 


#! /usrA)in/env python 

# QSL_Simons.py 

# Version 1.0 

# Copyright Jan-Ake Larsson <jan-ake.larsson01iu.se> 2015 
from Quantum_Siiiiulation_Logic in^xut * 

def Uf (state, s,peiiiO : 

n=state.n/2 # At this point, we have 2n elementary systems 
def Us(state,s): # For each s, construct one common U_s 
prev_i^one 
for i in xrange(n): 
if s[n-l-i]="l": 
if prev_i!^one: 

CN(state, (prev_i+2’tn,i)) 

CN(state, (i+2*n, i)) 
prev_i=i 
else: 

CN(state, (i+2*n, i)) 

state.extend(QSL_system(0,n)) # Add ancilla with ccxnputational bits 0 
Us(state,s) # U_s acts on input register and ancilla 

Peimute(state,perm,xrange(n)) # Peimutation on the ancilla 
for i in xrange(n): # Copy result from ancilla to output register 

CN(state, (i,i+n)) 

InversePermute(state,perm,xrange(n)) # Inverse peimutation on the ancilla 
Us(state,s) # Reset ancilla coipitational bits to 0 

state.remove(n) # Remove zeroed ancilla 

### Choose a bit mask s, and a randan peimutation (can be explicit and 
### non-random if desired, see below). A long s will work, but it will not be 
### possible to print all the function values 
#s="00000" 

#s="01000" 

s="01001" 

#s="lllll" 

#s="".join( [random.choice(("0","l")) for i in range(lOO)]) # 100-bit input 
#s="".join( [random.choice(("0","l")) for i in range(500)]) # 500-bit input 
#s="".join( [random.choice(("0","l")) for i in range(lOOO)]) # 1000-bit input 

n=len(s) 

### A random peimutation 
peinFRandPeim(n) 

### An e^qjlicit peimutation can also be used, for exanqDle mimicking a 
### 4-Toffoli as follows here. 

# pem=[0,1,2,3,4,5,6,15,8,9,10,11,12,13,14,7] 

### Note that the efficiency of Simon’s algorithm refers to uses of the 
### oracle, not the oracle storage or internals, here the peimutation. The 
### in^xiitant feature is that the oracle is used only linearly many times. 

print n,"bit input" 

### For small n, the function values can be canpared 
if rKlO: 

print "All function values" 
for i in xrange(2*’tn): 

state=QSL_syst€sn(i ,n)-H3SL_system(0,n) 

Uf (state, s ,perm) 

outcane=measureZ(state) # All bit indices 
print "f (“/os)=%s"“/o(outcome[0:n] ,outcome[n:2*n]) 
print 


print "Simon’s algorithm:" 

outcomes= □ 

basi^G 

while len(basis) < n-1: 

### Simon’s subroutine 
state=QSL_system(0,2=tn) 
for i in xrange(n,2*n): 

H(state,i) 

Uf (state, s ,pernO 
for i in xrange(n,2*n): 

H(state,i) 

y=mieasureZ(state,xrange(n,2=+n)) 

### End Simon’s subroutine 
outcomes. append(y) 

# remove linear span of outcomes from y 
for V in basis: 
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lp=v.index("l") # returns first index of a set bit 
if y[M="l": 
y+=V 

# Add y to the basis if y is still nonzero 
if y!="0"’tTi: 

# Remove the con^xment along y from previous basis vectors 
lp=y.index("l") # returns first index of a set bit 
for j,v in enumerateCbasis): 
if vD<]="l": 
basis [j] 
basis. append (y) 
if rKlO: 

print "Outcomes:", outcomes 
if rKlO: 

print "Outcone basis:", basis 
### s-star is orthogonal to the set of outcomes 
sstar="l"=+n 
for V in basis: 
if v.dot(sstar): 
sstar<-=v 

### Test the proposed s-star: is f(0)=^(s-star)? 
state=QSL_syst€m(0,2=+n) 

Uf (state, s ,penn) 

outccmeOmieasureZ(state,xrange(0,n)) 
state=QSL_syst€m(sstar,n)-H5SL_system(0,n) 

Uf (state, s ,penn) 

outcomesstar=measureZ(state,xrange(0,n)) 
if outccmeOl^utcomesstar: 


sstar="0"=fii 

print "After", len(outcomes),"oracle calls, computed s:",sstar 


print "Modified BH algorithm,", n, "oracle calls: 

outcomes [] 

for k in xrange(n): 

### Modified BH subroutine 

St at e=QSL_sy stem(0, n) -K5SL_sy stem( 1« (n- 1-k), n) 

for i in xrange(2*n): 

H(state,i) 

Uf (state, s ,perm) 
for i in xrange(n,2*n): 

H(state,i) 

outcomes.append(measureZ(state,xrange(n,2=tn))) 
### End modified BH subroutine 
if rKlO: 

print "Outcomes:", outcomes 
### s is orthogonal to the set of outcomes 
sstar="l"=+n 
for V in outcomes: 
if v.dot(sstar): 
sstar+=v 

print "Ccmpited s:",sstar 





